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We study the dynamics of ordering of the nonconserved and conserved Heisenberg magnet. The 
dynamics consists of two parts — an irreversible dissipation into a heat bath and a reversible preces- 
sion induced by a torque due to the local molecular field. For quenches both to T = and T = T C , 
we show that the torque is irrelevant when the dynamics is nonconserved but relevant when the 
, dynamics is conserved and is governed by a new nontrivial fixed point. 
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I. INTRODUCTION 



Most interacting systems like magnets and binary fluids when quenched from the disordered phase to zero temper- 
ature, evolve extremely slowly to establish order, primarily because of the slow annealing of the interfaces (defects) 
Ch ' separating the competing domains. It is seen that at late times, the system organizes itself into a self similar spatial 
distribution of domains characterised by a single diverging length scale which typically grows algebraically in time 
i ^i ' L(t) ~ t x l z . The equal-time order parameter correlation function C{r, t) = (0(r, t) ■ 0(0, t)) is a measure of the spatial 
distribution of the domains, and at late times is found to behave as f(r/L(t)), where L(t) is the distance between 
defects. The autocorrelation function, C(0,ti = 0,^) = (0(0,0) • 0(0, i2)), is a measure of the memory of the initial 
£^ , configurations, and decays at late times as L^t^) . The independent scaling exponents z and A and the scaling 
ly-j function f(x) characterise the dynamical universality classes at the zero temperature fixed point jjj. 

There has been a trend in recent years to compare the theories of phase ordering dynamics with numerical simula- 
tions of Langevin equations. Comparison with experimental systems, such as magnets, binary fluids or binary alloys 
have to take into account the various 'real' features that might be relevant to its late time dynamics. For instance, 
On : theories of binary fluids have to include effects of hydrodynamics, while those of binary alloys have to incorporate 
elastic and hydrodynamic effects. In the same vein, any comparison with the dynamics in real magnets has to include 
the effects of the torque induced by the local molecular field . 

s; 

^ . II. PHASE ORDERING DYNAMICS : QUENCHES TO T = 

g : 

The spins a (a = 1, 2, 3) in a Heisenberg ferromagnet in three dimensions experience a torque from the joint action 
of the external field (if present) and the local molecular field. In response the spins precess with a Larmour frequency 
^ 1 ilr, about the total magnetic field. Coupling to various faster degrees of freedom like phonons, electrons and magnons 
causes dissipation and an eventual relaxation towards equilibrium. 
The equations of motion in dimensionless variables is given by Q , 



^ = -HVf(v 2 + $- {$■$) 0) +g (0x V 2 0) • (1) 

The exponent /1 in the above equation takes the value when the magnetisation is not conserved (NCOP) and 2 
when it is conserved (COP). The dimensionless parameter g ~ fi_L/r is the ratio of the precession frequency to the 
relaxation rate, which is in the range g ~ 10~ 3 — 10 for typical ferromagnets. 
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A. Langevin Simulation 



We prepare the system in the paramagnetic phase and quench to zero temperature. We study the time evolution of 
the spin configurations as they evolve according to Eq. (0). We calculate the spatially averaged equal time correlator 
C(r, t) and the autocorrelator C(0,t — 0, £), both averaged over the random initial conditions. We compute the 
scaling exponents z and A and the scaling function f(r/L(t)) as defined above. L(t), a measure of the distance between 
defects, is extracted both from C(r, t) and from the scaling form of the energy density, £ = y J dr ( ( V0(r, i)) 2 ). 

We discretize Eq. ([!]) on a simple cubic lattice of size N ranging from 40 3 to 60 3 , adopting an Euler scheme for 
the derivatives fl(i[| with periodic boundary conditions. The space and time intervals have been chosen to be Ax = 3, 
At = 0.01 (NCOP) and Ax = 2.5, At = 0.20 (COP). Further details of the simulation may be found in 

1. Non-Conserved Case 

We will show that the torque is irrelevant in d = 3, which means that the asymptotic values of z, A and the scaling 
form f(x) remain unchanged and independent of g. The correlation functions have been computed for g ranging from 
to 5. The dynamical exponent is found to be unchanged from z = 2 (within statistical errors) and independent of 
g Q. Likewise the autocorrelation exponent is unchanged from its g — value and equal to A « 1.52 independent of 
g. Note that the numerical determination of A is subject to large errors and so we have to go to very late times 
and hence large system sizes. 
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Fig. 1. Scaling plot of C(r,t) = f(x) versus x = r/L(t) for g — 0(o), 0.5(A), 2(*). The solid line through the data 
points is a fit to the BPT function ||. Inset compares C(r, t) for g = l(o) with g = 5(D) at late times. The solid line 
through these data points comes from a preasymptotic analysis (see below) of C (r, t) g . 



The scaling function f(x) is also unchanged and independent of g. This is clearly apparent for g = 0, 0.5, 1, 2 (Fig. 
1). However at larger values of g, the correlation function crosses zero at large r, dips through a minimum, and then 
asymptotically goes to zero (of course / C(r, t) > 0). It would appear that C(r/L(t)) for g — 5 is qualitatively different 
from the scaling function for g — 0. However we note that the dip decreases very slowly with increasing time suggesting 
it might disappear at late times. Patience confirms this for intermediate values of g (between 2 and 5). When g = 5, 
finite size effects prevent the system from exploring its true asymptotic regime, when the order parameter field has 
totally relaxed with respect to defect cores. Pre-asymptotic configurations typically consist of spin wave excitations 
interspersed between slowly moving defects separated by a distance L(t) ^> £, the size of the defect core. Decomposing 
(j> into a singular (defect) part (p s i ng and a smooth (spin wave) part <j) sm gives the preasymptotic correlation function 
three contributions - C S i ng = (0sm 9 (O,t) • <j) s i ng (v,t)) (defect contributions), C sm = (0 sm (O,i) • (f) sm (r,t)) (spin wave 
contributions) and C SC at = (</>sm(0, t) • (f> S ing(r,t)) (scattering of spin waves off slowly moving defects). The defect 
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contribution C s i ng takes the form given in |5| , while the spin wave contribution C sm can be evaluated perturbatively 
[PJ. Even without including the scattering between spin waves and defects, we find that the total correlation function 
has a dip when g ^ as shown by the solid line in the inset of Fig. 1. The long lifetime of the dip arises because of 
the slow relaxation of the spin waves as they scatter off the slow moving defects. 



2. Conserved Case 



When the dynamics is conserved, we show that the the torque g drives the system to a new fixed point characterised 
by a different value of z, A and the scaling function f(x). Moreover these quantities are independent of the value of 
g as long as it is g > 0. This crossover is described by a crossover exponent and a crossover function. 

Figure 2 shows that the scaling function for g — is very different from those for g > and the g > scaling 
functions do not depend on the value of g. The z exponent crosses over from z = 4 (its value at g = 0) when t < t c (g) 
to z — 2 when t > t c (g) where t c (g) is the crossover time which decreases with increasing g (Fig. 2(inset)). Likewise 
the autocorrelation exponent A crosses over from A rs 2.2 when g = to A w 5.15 when g > ||. 
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Fig. 2. Scaling plot of C(r, t) comparing g = 0(o) with g > (g = 0.1(*), 0.3(A) 0.5(D). Inset shows log-log plot of L(t) 
versus t. At g — (o), z — 4 (line of slope 0.25 drawn for comparison). At other values of g (<? = 0.1(*), 0.3(D), 0.5(A)), 
z crosses over from 4 to 2 (line of slope 0.5 drawn for comparison). 



The crossover may be understood from a simple scaling argument On restoring appropriate dimensions, the 
dynamical equation Eq. (jl|) can be rewritten as a continuity equation, d<p(r, t)/dt — — V ■ j where the "spin current" 

y Q = -r(v^ + " eQ , 7 ^v^) . (2) 

Using a dimensional analysis where we replace j a by the 'velocity' dL/dt, we find 

§ = r£ + n^, ,3, 

where the parameters Mq, a and T _1 are the equilibrium magnetisation, surface tension and spin mobility respectively. 
Beyond a crossover time t c (g) ~ (Y/Mq£1) 2 ~ l/g 2 , simple dimension counting shows that the dynamics crosses 
over from z = 4 to z = 2 in conformity with our numerical simulation. Our numerics supports a scaling form 
L(t,g) = t 1 / 4 s(tg't') with w 1.7, valid for all g @. 
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B. Numerical Tests of the Mazenko Closure Scheme 



There is a class of approximate theories based on the gaussian closure approximation that is supposed to 
provide fairly successful predictions for the forms of the correlation functions. We will show that the gaussian closure 
approximation works fairly well for the nonconserved Heisenberg model with g > and leads to conclusions similar 
to the ones described in the previous section. In contrast, the use of the gaussian closure in conserved models leads 
to sharp inconsistencies. 

The gaussian closure method consists of trading the order parameter </>(r, t) which is singular at defect sites, for an 
everywhere smooth field m(r, t), defined by a nonlinear transformation, <j){v,t) = a (m(r, f)) . An appropriate choice 
for this nonlinear function a is an equilibrium defect profile, 

= j|^«' (4) 

where g(Q) = and 5(00) = 1. With this choice, m has the interpretation of being the distance away from a 
defect core. Correlation functions are calculated making the single assumption that each component of m(r,t) is an 
independent gaussian field with zero mean at all times ||[l| . 

To check whether this is a good approximation we compute <fi numerically and then determine m(r, t) by inverting Eq. 
(f|). We then calculate both the single point probability distribution P(m(r, i)) and the joint probability distribution 
in the scaling regime and compare with the gaussian assumption of Mazenko. These computations require a lot of 
averaging over initial configurations to obtain good statistics. 



1. Nonconserved Case 



Figure 3 shows the scaling plots of P(mi(r, t)) (mi is a component of m) at g — and g = 1. Since m has 
dimensions of length, the appropriate scaling variable is mi/L(t) The scaled distribution P(m\) is seen to be 
the same for g = and g = 1, suggesting that it is unchanged on addition of the torque. We also find that it is 
independent of g. 

The scaled distribution P(mi) closely resembles a gaussian at late times (this is confirmed by a more detailed 
analysis [Q]). The calculated joint probability distribution also agrees reasonably well with the Mazenko assumption 
H . These findings are consistent with a similar analysis done on a scalar order parameter . 




Fig. 3. Scaling plot of the un-normalised P(mi) for g = 0(t = 5000(o), t = 10000(D)) and g = l(t = 5000(A), t = 
10000(*)). 
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2. Conserved Case 



In the conserved case, the single point probability density P(mi(r, t)) scales at late times, but the scaled distribution 
is distinctly flatter than a gaussian (Fig. 4). Apart from this strong deviation we find that the gaussian assumption 
leads to an internal inconsistency, using a criterion initially developed by Yeung et. al. Q for a conserved scalar order 
parameter. 

We numerically evaluate the spectral density, the fourier transform of j(r,t) = (m(r + x,£) • m(x, t))/(m 2 (r, t)). We 
observe (^] that the spectral density, which should be a strictly positive function of its arguments, becomes negative 
for k/k m < 0.5 (j(k,t) is peaked at k m ) and in the range 1.5 < k/k m < 3.0 ! This demonstration highlights the 
intrinsic inconsistency of the gaussian approach for conserved vector order parameters. 




mi/L(t) 



Fig. 4. Scaling plot of the un-normalised P(mi) for g = at different times t = 900(o), t = 3600(D), t = 6300(A). 



C. Multiscaling Analysis 

We have just seen that the Mazenko assumption of a gaussian probability distribution fails sharply when the 
dynamics conserves the total magnetisation. This seems to lead to an interesting technical problem. Coniglio and 
Zannetti || had explicitly shown that the late time dynamics of a conserved n-component spin model in the limit 
ri — > oo reveals an infinity of length scales leading to a more complicated multiscaling form for the structure factor 
(fourier transform of C(r, t)) S(k,t) ~ L(t) p( ~ k / km ^ d . The length scale L(t) ~ i 1 / 4 and the position of the maximum 
of the structure factor A:" 1 ~ (i/lni) 1 / 4 , and the scale dependent exponent p(x) = 1 — (1 — x 2 ) 2 . They went on to 
speculate that this multiscaling behaviour might be a generic feature of conserved order parameters ||] . 

Subsequent numerical work showed that for the scalar (n = 1) ||, the XY (n = 2) [ ^o|JTT| and n = 3,4 |ll]] 
models, the structure factor obeyed the conventional 'single length' scaling form. This suggested the possibility that 
multiscaling was nongeneric and an 'aberration' of the n — > oo model, though there was no general proof. 

Bray and Humayun (BH) [|l2| provided such a 'proof. Their analysis was built on the validity of the Mazenko 
gaussian assumption for all values of n. Performing al/n expansion within the Mazenko framework, they showed 
that the asymptotic 5(k, t) exhibited multiscaling only when n — oo. 

Having just demonstrated that the Mazenko gaussian assumption is inconsistent for the conserved Heisenberg model 
(n = 3), how do we understand the BH proof []l 2f| ? Are the conclusions arrived by BH incorrect ? 

We now show that the structure factor 5(k, t) of the conserved n — 3 model does NOT obey multiscaling, consistent 
with the main conclusion of BH. This is done using the method described in |Io[ . 

The method demands a very accurate determination of S*(k, t). Since the numerical evaluation of S(k,t) 
is subject to large errors (especially at small k), we fit a function Cf(r,t) to the computed C(r,t) and then 
calculate the fourier transform Sf(k, t). The fitting function for Cf(r,t) has been taken as sm(r / L)/(r / L)(l + 
a(r I L) 2 ) exp(— 6(r / L) 2 ) which is similar to the analytic form given in |13j| . Note that only b and L are independent 
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fitting parameters, a is determined by the condition Sf(k = 0, t) = 0. We now plot Sf(k,t) versus t at fixed values 
of x = k/k m (Fig. 5). The resulting straight lines labeled by different values of x, all show a constant slope of 
approximately 3/4 (Fig. 5). Using the proposed multiscaling form, a plot of p(x) versus x (inset Fig. 5) shows that 
p(x) is clustered around 1. The small spread of p(x) around 1 indicates that we have not quite reached the asymptotic 
regime, and it is likely that the late time p(x) — > 1, in agreement with conventional scaling. In addition, note that 
the form of p(x) is qualitatively different from the downward curving p(x) predicted by Coniglio and Zannetti ||. We 
conclude then that the correlation function C (r, t) for the n = 3 conserved model does not obey multiscaling. 
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Fig. 5. Plot of S(k, t) vs t for x = 0.4(circle), x — 0.6(pentagon), x — 0.8(square), x = l.O(triangle), x — l.l(fillcd 
circle), x = 1.3(filled pentagon), x = 1.5(filled square), x = 1.7(filled triangle). A straight line of slope 3/4 is given 
for comparison. Inset figure shows plot of p(x) with x. Note that the errors in p(x) increases as | x — 1 | increases 
because of the smallness of 5(k, t) near its wings. 

Thus we are forced to conclude that since the Mazenko assumption is invalid for conserved vector order parameters 
(except when n = oo), an analytical proof of the absence of multiscaling of the asymptotic structure factor is still 
lacking. 



III. PHASE ORDERING DYNAMICS : QUENCHES TO T = T c 

We have also studied the asymptotic dynamics following a quench to the critical point || . We ask for whether the 
torque g is relevant at the Wilson-Fisher fixed point of the pure 5 = Heisenberg model, both when the magnetisation 
is conserved and nonconserved. 

Power counting reveals that g is irrelevant for d > 2 when the order parameter is nonconserved. Thus the values of 
z and A are unchanged from the Wilson-Fisher value. For the conserved model, on the other hand, power counting 
shows that g is relevant for d < 6, giving rise to a new fixed point. The dynamical exponent z changes from the WF 
value of 4 — rj to z — A — e/2 + C(e 2 ), where e = 4 — d and e = 6 — d p4| . The value of the autocorrelation exponent A 
is however solely determined by the conservation law; it is always fixed at the spatial dimension d. This is confirmed 
by a perturbative analysis to all orders in e || . 
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